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In a system with long-ranged correlations, the behavior of correlation 
functions is sensitive to the presence of a boundary. We show that surface 
deformations strongly modify this behavior as compared to a flat surface. 
The modified near surface correlations can be measured by scattering probes. 
To determine these correlations, we develop a perturbative calculation in the 
deformations in height from a flat surface. Detailed results are given for a 
regularly patterned surface, as well as for a self-affinely rough surface with 
roughness exponent Q. By combining this perturbative calculation in height 
deformations with the fleld-theoretic renormalization group approach, we also 
estimate the values of critical exponents governing the behavior of the decay 
of correlation functions near a self-affinely rough surface. We find that for the 
interacting theory, a large enough C, can lead to novel surface critical behavior. 
We also provide scaling relations between roughness induced critical exponents 
for thermodynamic surface quantities. 
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I. INTRODUCTION 



In a material with long-ranged correlations, such as a liquid crystal or a superfiuid, any 
local perturbation has influence over large distances. As a result, local properties, such as 
magnetization density, as well as correlation functions are modified on approaching a surface. 
Critical behavior near surfaces or defects, which is quite different from the bulk, has been 
extensively studied by means of the field-theoretic renormalizat ion-group approach HTj-Q. 
In this case, the local order parameter $ is perturbed near the surface up to a distance 
set by the diverging bulk correlation length C, ~ \T — Td'"^ , where Tc is the bulk critical 
temperature. Theoretical predictions for surface criticality have been tested experimentally 
and in simulations ||TO| , pJ]| . In particular, the grazing incidence of x-rays and neutrons 
^ has become a standard tool for probing critical behavior near surfaces and interfaces 
0-^. For instance, the decay of the two-point correlation function has been measured close 
to the surface of a FesAl crystal near its continuous order-disorder transition by the method 
of grazing incidence of x-rays f^. The phenomenon of critical adsorption near columnar 
defects has apparently been observed by small angle scattering of light in a NH4Br 
crystal near a continuous structural phase transition [[12 1. 

Most theoretical investigations so far have been restricted to flat surfaces. This is justified 
to a certain degree, since microscopic deviations from this idealized picture such as terraces 
of monoatomic height do not change the universal surface critical behavior [0,0 . However, 
for deviations on mesoscopic length scales, new phenomena are expected. Such deviations 
can be divided into two classes: 

(i) Advanced experimental methods of nanoscience such as x-ray ||15[, guided growth 
||T6| , and nanosphere lithography |T^, allow one to endow surfaces with specific, regular 
geometrical patterns down to the nanometer scale. These structures hold much promise for 
applications towards nanochips \TE\ or optoelectronic devices []T^. The surface modulations 
also offer a wide range of possible applications in fluid environments. For instance, at 
temperatures between the wetting temperature of the corresponding planar substrate and 
the critical temperature Tc of the bulk fluid, one can manipulate the adsorption properties of 
the fluid on the substrate by endowing the surface with periodic patterns of various shapes 

(ii) Surfaces or interfaces can be naturally rough, e.g., due to growth, fracture, or erosion. 
One possibility is that the substrate has a fractal surface, so that the surface area S grows 
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as a power of the projected area, i.e., S ~ L'^f where L is a characteristic length and 
df is the fractal dimension of the surface. Recently, the scaling behavior of correlation 
functions in a critical system in two dimensions near the fractal boundary of a random walk, 
for which df = 4/3, has been studied by methods of quantum gravity [^] and conformal 
invariance [Q. Another possibility is that the substrate has a self-ajjine surface, for which 
the surface area is proportional to the projected area. In this case the height fluctuations 
are characterized by a roughness exponent ( with < C < 1? so that (Sh)'^ ~ L^'' where 
6h is a typical height fluctuation over a distance L. Self-affine scaling is predicted by many 



numerical and analytical models of surface growth [p4|,p5| , and is also observed in a number 



of experiments |26|. A liquid- vapor interface, which exhibits rippled configurations due to 



the occurrence of capillary waves, is another realization of a self-affine rough surface ^Tj . An 
example where such an interface confines a critical system is given by the interface between 
liquid ^He near the normalfiuid-superfiuid transition and its noncritical vapor, which occurs 
in a recently used experimental setup in which the Casimir force in a critical system is 
measured ||2^ (see also Ref. p9|). 

In a previous Letter [Q, we showed that the shape of the surface has a distinct infiuence 
on the properties of an adjacent medium with long-range correlations. Here we demonstrate 
this in more detail for two-point correlation functions near a critical point of the medium, 
for both cases (i) and (ii) outlined above. Apart from Appendix 0, we choose the Dirichlet 
boundary condition $ = at the surface, which represents the so-called ordinary surface 
universality class in case of a fiat surface, and is usually appropriate for magnets, binary 
alloys near a continuous order-disorder transition, and '^He near the normalfiuid-superfiuid 
transition In Ref. |^T|, the infiuence of surface roughness on the fiuctuation properties 
of wetting films, and on the demagnetizing factor of a thin magnetic film, have been studied. 

In order to study the effects of the surface shape, we develop a perturbative expansion of 
two-point correlation functions in the deformations of the height profile. The method is the 
path integral approach used previously to calculate free energies |3^, and in the context of 
the dynamic |Q and static Casimir effect. Initially for a Gaussian field, the calculations 
are carried out to second order in the deformations. The first order results can also be 
derived by means of the stress tensor in conjunction with a novel type of short distance 
expansion (see Appendix |B|), and hold quite generally for any critical system bounded by 
a surface with either (a) Dirichlet boundary conditions $ = 0, or (b) boundary conditions 
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that break the symmetry of the order parameter near the surface. In the latter case, the 
leading singular behavior can be obtained by setting $ = oo at the surface, corresponding 
to the extraordinary or normal surface universality class, describing critical adsorption of 
a binary liquid mixture on the surface of a substrate or the interface between the critical 
liquid and its noncritical vapor The second order results are particularly useful for 

cases in which the first order contributions vanish (see below). 

The diffuse scattering of x-rays and neutrons at grazing incidence due to the modified 
correlations appears in addition to what would be observed if the surface was separating 



two homogeneous media |^5|]. The modified correlations may thus provide an additional 
and indirect means of characterizing the surface profile. This may be of value when other 
techniques are not possible, as in the case of the interior surface of a glass, or an internal 
crack, whereas scattering from a critical fluid or binary alloy coating the surface may be 
feasible. Already at the first order, the two-point correlation functions track the profile from 
the substrate, with a modulation that decreases with the distance of the two points from 
the surface. This leads to explicit predictions for the structure factor, as a function of the 
lateral wave vector transfer, for a modulated surface. 

For self-affinely rough surfaces, second order calculations are necessary, as the first order 
results vanish on average. In this context, the surface roughness is an example of quenched 
randomness. For a massless Gaussian field, we find the expected result that self-affine 
roughness leads to subleading corrections to the decay of two-point correlation functions, 
which at a scale r are smaller by a factor of r"^*^^~^) than the leading contribution coming 
from a fiat surface. Typical critical systems, however, are described by a non-Gaussian 
(interacting) field theory. In this case, the correlations are calculated perturbatively in a 
double expansion in the deformations and in the strength of the interaction, and the results 
interpreted with the aid of the renormalization group (RG) in 4 — e dimensions. We find that 
the subleading corrections now fall off with a slower power as compared to the Gaussian case 
and, surprisingly, for a sufficiently large ( even dominate, giving rise to novel surface critical 
behavior. However, for the XY model in two dimensions, below the Kosterlitz-Thouless 
temperature, we again find that the surface correlations fall off with the simple relative 
factor of r~^*^^~^) as compared to a flat surface (line). 

The results for correlation functions can also be related to thermodynamic quantities. 
To this end, we introduce distinct flelds hi, and hg in the bulk and close to the surface. 
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respectively, and propose a scaling ansatz for the leading singular part of the surface free 
energy per projected area By taking suitable derivatives of /J'^™^-' with respect to 

/ife and hg, we then obtain scaling relations for a variety of critical exponents related to 
thermodynamic surface quantities. 

The rest of the paper is organized as follows. In Sec.|I| we introduce the geometry, 
and develop the formalism for the perturbative calculation of correlation functions for a 
free (Gaussian) field theory. In Sees. |T| and we then consider a regularly patterned 
surface and a self-affinely rough surface in more detail. In Sec.|V] we combine the previous 
results with the RG, and obtain results for surface critical exponents. In Sec.jV^ we consider 
the XY model. Finally, in Sec. |VII| , we draw our conclusions and outline some possible 
extensions of our work; in particular, we relate our previous results for correlation functions 
to thermodynamic surface quantities via scaling relations. Some technical details are left for 
Appendices ^ - |C[ In Appendix for instance, we introduce a new type of short-distance 
expansion for the stress tensor. 

II. GEOMETRY AND FREE FIELD THEORY 

We consider a manifold Vl with the shape of a deformed surface. Each point on the 
manifold is represented by a vector X{y) = [X'^(y); ^ = 1, . . . ,d]; a D dimensional manifold 
Q embedded in d dimensional space is parametrized by y = (?/i, . . . , t//?). In the absence of 
overhangs and inlets, the surface profile can be described by a single-valued height function 
hly), where y spans a D = d — 1 dimensional base plane (see Fig. |l]). The parametrization 
of the surface is thus X(y) = [y^h^y)]. Position vectors r are decomposed according to 
r = (ry, z), where ry comprises the D = d — 1 components parallel to the surface, and z is 
the distance from the base plane. The vertical distance of r from the surface is given by 
S = z — h{r\\) (see Fig.|l]). We denote d dimensional vectors with underlined letters, and D 
dimensional vectors with boldface letters. 

Fluctuations in the critical system located above the surface will be described by an 
n-component order parameter field $(r) = [^i{r), . . . ,^n{ll)]- We consider the statistical 
Boltzmann weight e~^^ with standard Hamiltonian [||,|| 

= j^d'r |1(V$)^ + + i ' (2-1) 
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FIG. 1. Position vectors r = (ry, z) and r' = (rj|, z') of the two-point correlation function in the 
critical system located above and bounded by a deformed surface. The surface profile is described 
by the height function /i(x), and the vertical distances of r and r' from the surface are given by 
5 = z — /i(r||) and 6' = z' — /i(rj|), respectively. 

where tq ^ T — Tc and uq is the strength of the interaction. In this Section, we study 
the Gaussian theory, for which mq = 0. The volume V consists of the space available to the 
critical system. The above expression must be supplemented by a boundary condition on 
the surface. We choose the Dirichlet boundary condition $ = 0, representing the ordinary 
surface universality class. In this case, for n = 1 the order parameter $ can represent the 
magnetization in a uniaxial ferromagnet or the deviation of the composition in a binary 
alloy from the critical composition, for n = 2 the magnetization in a XY-magnet or the 
superfluid order parameter of ^He near the normalfiuid-superfiuid transition, and for n = ?> 
the magnetization in a Heisenberg ferromagnet 

The Gaussian two-point correlation function (or propagator) 

($,(r)$,(r')) = G(r; r') , = , (2.2) 

where the brackets ( ) denote the thermal average according to Eq. ( p.l| ) with = 0, can 



be calculated using functional integral methods p^j33[| . The details of this calculation are 
left to Appendix 0. The result is 

G(r;r') = G5(r;r') - j d^x j d^y (2.3) 
X G,[r; X, /i(x)] M(x, y) G,[r'; y, h{y)] , 
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where 

7£> 



G.fcr') = (2-4) 

with p = IpI, is the Gaussian propagator in unbounded bulk, and the kernel M(x, y) is the 
inverse of the kernel (^^[x, /i(x); y, h{y)], i.e. 

I rf^i/M(x,y)G,[y,My);y',My')] = 5^(x-y'). (2.5) 

While the above results [with an appropriately modified bulk propagator in Eq. (p.4|) ] are 
generally valid, we focus on the behavior of the correlation functions at the bulk critical 
point, i.e., for T = Tc, where correlations are strongest |]36| . 

Equation ( |2.3| ) is difficult to evaluate in general. To proceed, we now consider the height 
profile h{x) as a small perturbation, and expand G{r;r^) in a series Go + Gi + G2 + ■ ■ ■ in 
powers of h up to second order, under the constraint that z and z' are kept fixed. The lowest 
order result, 

Gofcr') = Gb{r\\,z-4,z')-Gb{r\\,z;r[,-z') (2.6) 



(2^ 



e 



*P-(r||-r;) 



^^-^^^>go{p;z,z') (2.7) 



with [see Eq. (|]4| 



goip; z, z') = 1 [e-^'^-^'l - e-^^^+^'^l (2.8) 
2p L J 

corresponds to a fiat surface, and can be obtained by the method of images ^j^. The bulk 
correlation function Gb{r;,r^) decays as r~^'^~'^~^^^ for large separations r = |r — Zl'I, where 
the bulk critical exponent rj is given by 77 = in the Gaussian theory. In contrast, if both 
points remain close to the surface, Godisr') decays as r^*-'^^^+''ii-*, where r/y is a surface critical 
exponent given by r/y = 2 in the Gaussian theory |l| 
The first order result is given by 



with 



Giinr!) = j rf^x J(r||,x;2)/i(x)J(rj|,x;z'), (2-9) 



J(x,y;z) = - ! e'P<^-y^ e~P' . (2.10) 



2 J (27r) 
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Note that J(x, y; z — O"*") = |5'^(x — y), where (5^(x) is the delta function in D dimensions. 
Already the result at this order tracks the profile /i(x) of the surface. For example, for 
p = |r|| — rjil ^ oo with z and z' fixed, the above results for Go and Gi imply the behavior 
(see Appendix 0) 

G(r; r') ~ [1 - A{r) - A(r')] p-^'^-^+^'il) , (2.11) 

up to terms of order (h/z)^ and (h/z'Y. Thus, the leading power law is the same as for a 
fiat surface, but the amplitude is modulated by the surface deformations in the vicinity of 
ry and rj| by the dimensionless and universal amplitude 

A{r) = ^ I rf^a:^A(x-r||,^) , (2.12) 

where for the present Gaussian case, 7]\\ — rj = 2 and A(x — ry, z) = 2J(x, ry; z). Equations 
( p. lip and (|2.12|) are valid quite generally, and in particular also for the boundary condition 



representing critical adsorption of a binary liquid mixture (see Appendix 0). The explicit 
form of A(x, 2;), however, depends on the surface universality class considered. 
The second order result reads 

G'2(r:;r') = J d^'x J rf^y /i(x)/i(y) C(r, r'; x, y), (2.13) 



with 



C(r,r';x,y) = - 8 J(r||, x; z) J(r;|, y; /) K(x, y; ^ 0+), (2.14) 



K{x,y;z) = - [ -^^^e'P-("-y)pe-P" . (2.15) 



and 



2 J (27r)« 

In a scattering experiment with grazing incidence, one probes the lateral structure factor 
^(p, z; p', z') 0J^,^, which is defined by the Fourier transform 

G(r;r') = / / e'-''i S(p,.;p',/) . (2.16) 

Using the Fourier transform of the height profile 

My) = / (l^e^'^'Mk), (2.17) 

with /i(— k) = h(ky, we obtain an equivalent expansion S = Sq + Si + S2 + ■ ■ with 
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2p 



^-p\z-z'\ _ ^-p\z+z'\ 



(27r)^5(p + p') , 



Si = -e-P'e-P''' h{p + p') 
S, = -e-P^e-P''' ■ 



p - k| h{k) h{p + p' - k) 



(2.18) 
(2.19) 
(2.20) 



(27r)^ 

For a rough surface, the deviations in height from a planar surface have no upper bound. 
In this case, it is convenient to carry out the expansion in h{x) for fixed vertical distances 
6 = z — h{T\\) and 6' = z — h{r'^^), instead of for fixed z and z' (see Fig.|I|). This representation 
facilitates the perturbative analysis of the field theory described by Eq. (|2.1|) (see Sec. 0). 
Moreover, in view of probing correlation functions lateral to the substrate surface by grazing 
incidence scattering of x-rays and neutrons [3,|6l,p5| , this representation is natural, since in 
these experiments 6 and 6' show up as length scales which are set by the finite penetration 
depth of the x-rays. 

Writing G = Gq + Gi + Gu + . . . where the subscripts 0, I, II indicate the corresponding 
order in h{x.) under the constraint that 6 and S' are kept fixed, we find 



Go(r;r') = Gfe(r||, 5; r' , 5') - Gfc(r||, 5; r' , -5') 



Gi(r;r') = - [h{r\\) - h{r[)] ^ G,{r\\,6;r\,6 



(2.21) 



(2.22) 



+ 2 Jd^'x J(r||,x;(5) [h{r\^) + h{r\^) - 2h{^)] J(rj|,x;5') , 



G'ii(r;r') 



K(r||,rj|;|(5-(5'|) + ir(r||,rj|;(5 + (5') [h{r\\) - h{r\)Y (2.23) 
J d^'x J d^y J(r||,x;5)Mo(x,y) [/i(x) - %)] V(rj|, y; 5') 

j d^x K{vii,x; 6) [/i(r||) - h{x)f J(rj|, x; S') + (r ^ r') 



The first line in Eq. ( p.22|) is valid for 6' < 6, and Mo(x, y) in Eq. ( p.23|) is defined as in 
Eq. ( |2.5D but with /i(y) = 0. The kernels J and K are given by Eqs. (|2.1CI|) and ( p.l5[ ), 
respectively. The contribution Go iii Eq. ( ^.211 ) corresponds to the Gaussian propagator for 
a half-space bounded by a flat surface with Dirichlet boundary conditions, i.e., 

d^P ip-(r||-r;) 



Go(r;r') 



(2vr) 



D 



(2.24) 



with go from Eq. (E 
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III. MODULATED SURFACES 



We now apply the results of the previous Section to patterned surfaces. The simplest 
example is provided by an uniaxial sinusoidal modulation with wavelength A along, say, the 
X direction, i.e., 

h{x,Y) = acos{27Tx/X) . (3.1) 

The other D — 1 directions along the surface, denoted by Y, remain translationally invariant. 
The Fourier transform of this height profile is 

h{k) = (27r)^ 6^-\K) [6{K - f ) + + f )] , (3.2) 

where k is decomposed according to k = {kx-, K). 

The nontrivial orders of the expansion of (^(r; r') in h for fixed z and z' are given by 

x-y / I l\ ^ ^^x' f jp(rii-r') /r, ^\ 

Gi(r||,z;r||,2;) = --e / ,^ e " (3-3) 



2 J {2tt)d 



G,{r\\,z;rlz') = - e^^'^^ii-^liMp - (f , 0)| e"^^ e'^^' (3.4) 

ip-(r||-r' ) 



/A ^ P jp(r||-r(|) 



a iiiix' P 



-e ^ 



4 J (27r)« 



X |p- (^,0)|e-f^e-|P-(^'0)l^' + (rii^rj. 



For p = |r|| — r'l I —>■ oo, the leading power law in p is the same as for a fiat surface, but the 
amplitude is modulated by the shape of the surface in the vicinity of ry and rj|. In particular, 
the first order result in Eq. (|3.3| ) is consistent with Eqs. ( ^.11 ) and ( |2.12| ). For z, z' <^ a, X 



the correlations follow more or less the surface modulation. Interestingly, for z, z' ^ X, the 
correlations that are sensitive to the modulation, i.e., depend on A, decay exponentially in 
z/X. For instance, for z = z' and z/X ^ oo, one has Gi ~ e~^^ and G2 ~ e~~^. This 
exponential decay is due to the fact that the surface profile has a perfect periodic 

shape. In contrast, a local perturbation on the surface would result in a perturbation of the 
correlations which decays only as a power law with the distance from the surface. 
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The corresponding orders of the lateral structure factor are given by 

^i(p, z- p', z) = -^ e-^^ e-^'^' (27r)^ ^^-^(P + P') (3.5) 
X [%. + p;-f) + 5(p.+P^ + f)] , 

52(p, p', z') = -'^ e'^-^ e'^'^' (27r)^ ^^-^(P + P') (3.6) 
X {|p - (f , 0)1 [5{p^ + - f ) + 5(Px. + pD] 
+ Ip + (f , 0)1 [5(p. + + f) + '^(P- + P^)] } • 

These results indirectly characterize the surface in scattering experiments. For instance, the 
form of 5*1 implies that the incident wave vector component px is scattered to p'^ = ± ^ 
while the other components of p remain unchanged. The form of 5*2 implies that is 
scattered by 0, — ^ scattering experiment with grazing incidence, the length 
scale perpendicular to the surface is set by the depth b that the evanescent wave penetrates 
the sample, giving rise to diffuse scattering and thereby probing the critical correlations 
close to the surface 0. Since this diffuse scattering appears in addition to the contribution 
already present away from criticality ||3^, it can in principle be separated out by tuning 
the temperature deviation T — Tc. We assume that b is much larger than the height of the 
deformations. In this case, the above expansion in the deformations results in an expansion 
in powers oi h/b for the elastic scattering cross section, which allows one to distinguish 
the corresponding contributions via their intensities. 

IV. ROUGH SURFACES 

The second order results are particularly useful when dealing with rough surfaces, where 
the quench averaged first order corrections vanish. Within the description using a height 
function /?.(x), self-affine roughness is characterized by the behavior 

[/i(x)-My)]2 ~ |x-y|2f , Ix-yHoo, (4.1) 

where the overbar denotes averaging over self-affine realizations of the surface profile, and C, 
with < C < 1 is the roughness exponent. Without restriction of the generality we choose 
the coordinate system so that /i(x) = 0. In the limit of short distances |x — y| it is reasonable 
to assume that the surface is smooth. This can be modelled by the Fourier transform 
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While at large separations the above correlations grow as |x — yp, we have also introduced 
a cutoff length A to regulate the behavior of the surface at short distances, and an overall 
amplitude length uj. The length A characterizes the crossover from the analytic behavior for 
|x"~y| ^ A to the behavior in Eq. ( |4. 1| ) for |x — y| ^ A. Apart from its physical significance, 
the appearance of the finite crossover length A in Eq. ([4.2|) is also essential within the present 
theoretical approach (see Sec.^). 

A characteristic feature of self-afiine roughness is statistical translational invariance, since 
the right hand side of Eq. ( |4.2| ) depends on the distance |x — y| only. This implies that the 
averaged lateral structure factor S is proportional to 5^(p + p'), and depends on z, z', and 
p = IpI only. In order to maintain translational invariance, it is convenient to express the 
results for the correlation functions in terms of the local distance 5 = z — h{Y\\) from the 
surface rather than z (see Fig. |1]). The two-point correlation function must now vanish as 
5 or 5' go to zero. On averaging G{r]r^) over different surface profiles, the contribution 
Gi in Eq. (|2.22|) vanishes due to h{x.) = 0, and the contribution Gu in Eq. ( p.23| ) becomes 
translationally invariant with respect to the lateral components ry and rj|. We thus introduce 
the lateral Fourier transform 



}C{p,6) is the lateral Fourier transform of -ft'(x, y;(5) [/i(x) — h{y)]'^ and we have used the 
fact that the lateral Fourier transform of Mo(x, y) [/i(x) — /i(y)]^ appearing in the second 
line of Eq. ( |2.23|) after averaging is given by 4/C(p, 5 = 0). Using Eq. (^3), IC{p,6) can be 
expressed in terms of the convolution integral 




(4.3) 



where (72 (p; ^, can be read off from the right hand side of Eq. (p. 23 ), i.e.. 



92{p;S,6') = l[/C(p,|5-5'|) + /C(p,5 + 5') 

+ /C(p, 0) e-P(^+^') - JCip, 5) e-P^' - ]C{p, S') e'^^ . 



(4.4) 



k:{p,s) 



= UJ' 



t/(lp-kl,5) 



(4.5) 



where U{p, 6) is the lateral Fourier transform of K{'x, y; 6) |x — yp given by 
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U{p,5) = 
In terms of the coordinates r 



1 



p5^ 



D 



1 /I 



6 



-p5 



(4.6) 



2' 2 

(r||,5), the above results imply that the leading power law 
oo is the same as for a flat surface. The corresponding 



■II 



2(1-01 



as 



behavior of G{r; r') for p 

amplitude depends on the roughness, and is modified by a factor of [1 — k{uj/\) 
compared to a fiat surface, where k > is a number of order unity. The subleading correction 
of order /i^ decays with the separation p with an additional factor of p"'^^^^'^^ compared to 
the leading term [see Eqs. (|5.16|) and ( ^.17|) in Sec. for e = 0]. 

Note that (72 (p; vanishes for 5 = or 5' = as it should, according to the Dirichlet 
boundary condition at the surface. This would not be the case for z = or z' = if 
we carried out the expansion in h{x) with fixed z and z'. However, the realization of the 
Dirichlet boundary condition for the Gaussian propagator is essential for the perturbation 
theory of the field theory described by Eq. ( |2.1| ). Moreover, (72 (p; S') is an analytic function 
for small 6 or 6' due to the finite crossover length A in Eq. ([4.5|) , which would otherwise be 
ill defined for A = if 5 = and C < 1/2. 



V. INTERACTING THEORY 



In this Section we consider the asymptotic scaling behavior of the two-point correlation 
function near a self-affine rough surface for the n- vector model at the bulk critical point. 
By combining our previous results with the field-theoretic renormalization group (RG), we 
estimate the values of the corresponding critical exponents, using a double expansion in the 
surface deformations and in the deviation e = A — d of the space dimension d from the upper 
critical dimension. 

For the interacting field theory, governed by Eq. ( |2.1| ) with Mq 7^ 0, standard perturbation 
theory can be applied to get the correlation function near a surface of arbitrary but fixed 
shape, 

($,(r)<l>,(r')) = 6,,g{r,r:;uo), (5.1) 

with 

^(r,r';«o) = G(r;r') - ^ y G{r; R) G{R- R) G{R:,r:) + 0{ul) , (5.2) 

where the Gaussian propagator G{'r^r[) is given by Eq. (|2.3| ). We are interested in the 
behavior of ($i(zL)$j(r')) limit for which the distance between r and r' is much larger 
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than one or both of the vertical distances S and S' (see Fig.|l|). If S' is small, say, it is helpful 
to consider the so-called surface operator 

<l>^(rj|) = 9„<l>(r') , (5.3) 

where (9„ = [(7(rj|)]^-'^/^[95/ — V/i(rj|) ■ V] denotes the normal derivative at rjj on the surface, 
with the determinant gly) = 1 + [V/i(y)]^ of its induced metric [see Eq. ( AlUI )]. In this way 



one avoids to deal with the irrelevant length 6' from the outset. For correlations vertically 
away from the surface, i.e., ry = rj|, we are thus led to consider 

{MD'^fi^)) = s.,gAr\\,S;uo) . (5.4) 

The loop expansion of g±{r\\, 6; Uq) is obtained by taking the normal derivative at ry of the 
right hand side of Eq. (|5l^) . 

Up to now in this Section we have considered a surface with arbitrary but fixed shape. In 
particular, for a flat surface, the one-loop addition in uo can be regularized and renormalized 
by minimal subtraction of poles in e = 4 — d, leading to logarithmic contributions in the 
separation r = |r — r'|. This perturbative result can then be improved by RG, resulting 
in power laws in r with corresponding surface critical exponents For a self-affinely 

rough surface, the function Q± depends, of course, on the shape of this surface, i.e., on the 
height function h{x.). However, by averaging over different surface profiles, we expect that 
the average Q± depends only on gross features characterizing the surface configurations, and 
in particular becomes independent of ry due to translational invariance. In the following 
we restrict ourselves to surfaces which are rough on large distances, and to contributions 
to Q± up to second order in /i(x). According to Eq. ( [4.2| ) we conclude that in this case 
the amplitude uj and the crossover length A are the only remaining relevant length scales 
characterizing the different surface configurations. 

In the next step, the resulting average Q±{6;uj,X]Uq) has to be renormalized. For our 
perturbative calculations we use dimensional regularization and renormalization by minimal 



subtraction of poles in e = A — d . The reparametrizations 

Uq = IQtx"^ if ZuU (5.5) 

and 

$ = 4/'$,e„ (5.6) 
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of the bulk parameter uo and the bulk field $ in terms of their renormalized counterparts 
u and $ren are not affected by the presence of the surface. Here = 1 + 0{u) and 
Zij, = 1 + are the corresponding renormalization factors, and /i is the inverse length 

scale which determines the renormalization-group flow. Since all surfaces we average over 
are smooth on short distances, i.e., distances much smaller than the crossover length A, we 
expect that the surface operator is renormalized by the same renormalization factor Zi 
which would occur for a flat surface with Dirichlet boundary conditions. Thus, 

$^ = iZ^Z,y/'^i^, (5.7) 

with [§ 

Zi = 1 + — - - + Oiu^) . (5.8) 

Using the above reparametrizations the renormalized, i.e., pole-free, counterpart of Q± is 
given by 

G±,ren{^',^,^',%IJ') = Z^^^'^ Q±{6; U, X] Uq) ■ (5.9) 

This perturbative result can be improved using standard renormalization-group methods, 
by noting that Q± does not depend on fi. The asymptotic scaling behavior is governed by 
the infrared (long-distance) stable flxed point for which 

u = u* = + 0{e^) , (5.10) 

and G_L^ren assumcs the scaling form 

^,,,J5;u;,A;w,/i) ~ r^-^-^^"^) hi5/X;uj/X) (5.11) 

with the critical exponent r]± for a flat surface. The scaling function f± is universal, but 
depends on the particular way we have introduced the crossover length A in Eq. ( [4.2|) . Since 
all surfaces we average over are smooth on short distances, f^{0;uj/X) should be a flnite 
number (in the following we suppress the dependence of f± on uj/X for brevity). In the 
other limit 6/X ^ oo, the scaling function /_l(5/A) is expected to exhibit a power law which 
reflects the self-affine structure of the surface. 

We have conflrmed Eq. ( |5.11|) explicitly to one- loop order according to Eq. ( |5.2|) , using 
the expansion of ^(r, r') up to second order in h{x.) in Eqs. (|2.21|) - (|2.23|) , and averaging 
using Eq. (|4.2| ) . Figure |^ illustrates this double expansion in graphical form . 
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FIG. 2. (a) Representation of the full Gaussian propagator G{r;r') in Eq. (|2.3|) and its expansion 
up to second order in /i(x) according to Eqs. ( p.21 ) - ( |2.23| ). The number of ticks corresponds to 



the order in /i(x). (b) The second order in /i(x) contribution to the one-loop integral in Eq. (|5.2| ) 
decomposes into several parts. The dashed lines connecting the ticks indicate averaging over 
different surface profiles, using Eq. (|4.2|). The cross corresponds to the surface operator 



We indeed find that the 1/e poles generated by the surface operator in Eq. 
are removed by the renormahzation factor Zi in Eq. ( [5.8| ), which provides a test of our 
calculation, and for the reasoning leading to Eq. (|5.11|) . This calculation gives also the 
explicit form of the scaling function f± to first order in e. We confirm, in particular, that 
/_l(0) is a finite number, and that the logarithmic contributions of f±{S/X) for 6/X —>■ oo 
can be recast in the form of a power law, i.e., 

f±{5/\) a + (3{5/\y . (5.12) 

Whereas both amplitudes a and (3 depend on u/X, the universal exponent ip is independent 
of u/X and given by 

^. = -^.-(2-20 +0{e'). (5.13) 
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Perpendicular correlations are obtained when r moves into the bulk, while r' remains 
close to the surface, i.e., 6 ^ oo with 6' fixed (see Fig. 0). Equations (|5.4| ) and ( ^.11|) - (|5.13|) 
then imply that the correlations decay as 

1 



iMrJMr')) - ^jz^ + , (5.14) 

where the first term corresponds to a fiat surface with ri± = 1 — + C>{6^) [|T],3. The 

second term describes the effect of self-affine roughness, with an amplitude a depending on 
uj, X, and (, and the new universal exponent 

77-1-2 

V± = V±-i^ = (2-20 + 1 - 2^5 + 0{e') . (5.15) 

72-1-8 

Similarly, when both points remain close to the surface, i.e., p = \r\\ — rj|| — ^ oo with 6 and 
6' fixed, the correlations decay as 

1 



iHrJMr')) - + . (5.16) 

In this case the fiat surface is governed by rj\\= 2 — e + O(e^), while self-affine roughness 
gives 

77-1-9 

^11 = (2-20 + 2-4^6 +0{e'). (5.17) 

77-1-8 

The corrections due to roughness now decay with a slower power as compared to the 
Gaussian case. Indeed, for a sufficiently large roughness exponent these corrections 
can even dominate the result for the fiat surface. The borderline roughness exponent is 
CI = l- i^^ + ^'(e^) for perpendicular, and a different value of Cii" = 1 - | ^ s + O^e^) 
for parallel correlations. This is a surprising result from a naive point of view since, due to 
C < 1 , on larger and larger length scales a self-affine rough surface looks more and more like a 
fiat surface. Note that this effect becomes only visible beyond the Gaussian approximation, 
which corresponds to e = 0. By setting e = 1 in the above expressions, one obtains the 
corresponding estimates for ci = 3. 
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VI. TWO DIMENSIONAL XY MODEL 



To compare the results of the previous Section with a different interacting theory, we 
examine the correlations for a two dimensional XY model below the Kosterlitz-Thouless 
temperature |40]. The order parameter in this system is the spin variable s{r) = e'^^^-\ 
where 6{r) is the angle (phase) the spin makes with some reference axis. Even though the 
phase fluctuations are described by a Gaussian model, nontrivial spin-spin correlations are 
obtained. Below the Kosterlitz-Thouless temperature, the d = 2 dimensional XY model is 
well described by the spin-wave Hamiltonian (neglecting vortices) I^TJ 



f3n{e} = \k j dPr{V9)\ (6.1) 
where r = {x,y). Correlation functions decay as power laws in this system. For instance. 



the two-point correlation function in the unbounded plane is given by [11,12 

Gb(r;r') = {e'^(r.)e-^oir')^ ^ exp[^5(r;r') - '^g,ir;r) - i6;,(r';r')] (6.2) 

with 

Gb{r;r:) = (^(r)^(r')) = " ^ ln(r/a) , (6.3) 
where r = |r — r'| and a is some lattice cutoff. This implies 

G.(r;r') = (-)'' , (6.4) 



-a/ 

where t] = 1/{2'kK). 

If the plane is bounded by a free surface (line) at ?/ = 0, the correlation function (^(r; r') 
in the half-space ?/ > is given by similar expressions as in Eq. ( |6.2| ), where now ^(r, r') 
satisfies the Neumann boundary condition at the surface The final result 



G{x,y;x',y') 



(6.5) 



Ayy' 

implies the surface critical exponents r/y = 2r] and r]± = ^r], which fulfill the scaling relation 
2ri± — ri^\ = 1] familiar from the surface critical behavior of n- vector models []T],0|. 

In order to study whether the nontrivial roughness dependence of correlations obtained 
in the previous Section is also present here, we now consider a deformed surface (line) with 
the same boundary conditions as above. Similar steps as outlined in Appendix ^ lead to 
the result for the two-point correlation function 
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G{r',r!) = exp[r(r;r') — |r(r;r) — |r(r';r')] (6.6) 

with 

r(r;r') = Gbinr:) - j j dx'dnQb[T:.X{x)]M{x,x')dn'Q,[r^;X{x')] , (6.7) 

where 9„ denotes the normal derivative acting on X, and M.{x, x') is the functional inverse 
of dndn'Gb[X{x),X{x')]. As in Sec.^, we use the representation X{x) = [x,h{x)] in terms 
of the height profile h{x), and expand G{r;r^) up to second order in h. In particular, for a 
self-affinely rough surface, we find, using Eq. ( [4.1|) , that the surface correlations fall off with 
the simple relative factor of r^^*-^^''^ as compared to a flat surface (line) (compare SecfTV]). 
We attribute this to the Gaussian nature of the fluctuations in the phase angle 6{r), which 
are retained in the asymptotics of correlations for s{r). 

VII. CONCLUSION AND OUTLOOK 

We have developed a path-integral formulation for the study of correlation functions 
in a system which is confined by a deformed boundary. Our results are generic for any 
system with long-ranged correlations. Examples include systems with a broken continuous 
symmetry, such as the XY model below the Kosterlitz-Thouless temperature, or a nematic 
liquid crystal, where the correlations are generated by the corresponding massless Goldstone 
modes; or critical fluids or magnets described by the n-vector model at the bulk critical 
point, which has been mostly considered in this work. The surface deformations can consist 
of specificly designed, regular patterns, or represent a self-affinely rough surface. Some 
conclusions and possible extensions of this work are listed below. 

(i) Thermodynamic surface quantities: Thermodynamic quantities can be obtained from 
derivatives of the free energy with respect to magnetic fields. To discuss surface behavior, 
we introduce distinct fields and hg in the bulk and close to the surface, respectively. 
Assuming that our underlying assumption of the validity of an expansion in h{x.) holds, the 
results for the two-point correlation function are consistent with the following form for the 
scaling of the leading singular part of the surface free energy per projected area. 



j(smg) _ ^-^^+1 



9s{he\h,^y^) + r'^''^^9r{h^'\h,e') , (7.1 



where ^ ~ |T — is the correlation length that diverges at the critical point. The first 

term in square brackets corresponds to a fiat surface, with i/b and ys describing the relevance 
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of bulk and surface fields, respectively |T|,0]. The second term gives the effect of surface 
roughness, with ^"^(i-c) reflecting the average increase in area. 

By taking derivatives of Eq. (|7.1| ), one can derive scaling relations between various surface 
critical exponents, in complete analogy to the case of a fiat surface . In the following we 
focus on the contributions generated by the surface roughness, which according to Eq. 
appear in addition to the corresponding contributions for a flat surface. For example, the 
singular part of the surface magnetization^ —dfs^^^^^/dhs, can be written as mi + mi so that 
Xi = drhi/dhh and Xn = dmi/dhg represent the contributions to the local susceptibility 
and the layer susceptibility generated by the surface roughness, respectively. Similarly, we 
suppose that the singular part of the two-point correlation function near the surface can be 
written as G{r_]r^) + G{r;r^), and behaves for hi, = hg = a.s 

f lz:-z:'r^°'~^+^"^ r||(|r-r'|/0 , ^ = o, 

G(r;r') ~ <^ (7.2) 

where is the angle r^ — r' makes with the surface, and r_|_ vanishes for -i? ^ 0. Equations 
1|) and ( [7.2|) then imply the scaling relations between various critical exponents related 
to a rough surface shown in Table |. 

Equations (|5.15 ) and ( ^.17 ) for the n- vector model are consistent with the scaling relation 
for rj\\ and rj± shown in Table |. However, to regain the results in Eqs. ( |5.14|) - ( ^.17| ), 
we have to use a value of = 1 + . e + 0{e^) in Eq. ( [7.1|) , which is different from 



Us = 1 — -Tg£^ + C^(^^)- To motivate and justify this difference, we resort to an analogy 



2(n+8) 

c J- n(c'^\ Tr. mril-lirQl- 

n+8 

in which the rough surface is replaced with a collection of edges with a (possibly scale- 
dependent) distribution of opening angles. Already for a single edge, describing correlations 
requires a distinct and angle-dependent value of ije for the magnetic field close to the edge 
43| , |4^ . Similarly, results obtained recently for correlations in the vicinity of a fractal surface 
with fractal dimension dj [^,0 cannot be obtained using the value of ys for a flat surface 
[with ^"'^■f replacing ^^"^+1 in Eq. ( [7.1| ) and omitting the second term in square brackets] . 
Thus ys can be regarded as inherently related to self-affine geometry. Interestingly, however, 
ijs itself does not depend on the roughness exponent C, at least to order e. 
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TABLE I. Scaling relations between critical exponents relevant to a rough surface, as derived 
from Eqs. ( |7. 1| ) and ([7.21), in terms of the bulk critical exponents r/, i/, = A/z/, and the 
roughness exponent For each exponent in the left column there is a corresponding exponent 
for a flat surface [1,2] which would be denoted without tilde (compare with Table III in Ref. [2]). 



critical exponent 



conditions 



scaling relation 



^x, ^11 Eqs.(|l|), ill) 
Vs Eq. (P) 



Xl ~ \T 



-71 



Xll ~ \T\ 



nil [—' 



-711 



mi 



l/<5i 
/ijV^ii 



T = hfy = hg = 

T^O,hb = hs = 
T^O,hb = hs = 
T < 0, hb = hs = 

T = hs = 0,h^O 
T = hb = 0,hs^O 



2^±-^|| =?? + 2-2C 
= i((i-^ll+2-2C) 
71 = 1^(2 - rj±) 
711 = u{l - ) 
^1 = f(d-2 + ^|| + 2-2C) 

(^ii = i^Vs/A 



(ii) Higher orders of the perturbation theory: As the previous remark already indicates, 
higher order results in e are necessary in order to check the generality of our results for the 
n- vector model. For the contributions up to second order in h{x.) (as considered here), we 
expect a systematic expansion in powers of e, and one can calculate the 0(£^) contributions 
of, e.g., fj_ and ip in Eq. ( ^.12| ). All the information needed about the self-affinely rough 
surface is contained in Eq. ( [4.2|) . However, it is not clear how the perturbative calculation in 
/;,(x), for a self-affinely rough surface, can be generalized to higher orders than the second. 
Such an attempt would require, in addition to Eq. (^]2|), the knowledge of stochastic averages 
of three and more fields h{x.), which can also introduce new length scales. Regarding these 
obstacles, it would be desirable to complement our results with a nonperturbative approach, 
e.g., for the two-dimensional Ising model bounded by a self-affinely rough boundary. 

(iii) Multiscaling: For a random fractal boundary, it has been shown [^] that correlation 
functions exhibit multiscaling, which means that the average (over fractal realizations of the 
boundary with given fractal dimension dj) of their nth power does not scale in the same 
way as the rath power of their average. It would be interesting to see if similar behavior also 
applies to self-affine rough boundaries. 
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APPENDIX A: PATH INTEGRAL FORMULATION FOR CORRELATION 

FUNCTIONS 



In this Appendix we introduce a method to evaluate correlation functions of a fluctuating 
field subject to boundary conditions at surfaces of arbitrary shape. We consider a scalar 
field $ described by the Gaussian action 

'1 



<fr 



2^^^)^ - 2 



(Al) 



corresponding to Eq. (p.l|) with uq = 0. In order to study the behavior of correlation 
functions for cases where more than one boundary surface is present, we consider manifolds 
(objects) Qa with a = 1, . . . ,N. Each point on the manifold Qa is represented by a vector 
^a(y) = (A^a(y); = l, ■ ■ ■ ,d). Assuming the Dirichlet boundary condition $ = on the 
manifolds, a general correlation function with respect to the action ( |AT| ) can be written as 

1 r ^ 



where 



a=l Xc 



N 



'' a=l Xa 



(A3) 



Correlation functions of $ can then be deduced from the generating functional 



Z{J] 



exp 



d^rJ{r_)^{r_) 



(A4) 



which is normalized such that ^{0} = 1. 

Following Refs. [^,0, we now express for each manifold Vt^ the boundary condition 
enforcing functional (5[$(Xq)] in terms of an auxiliary field '^a{Xa) as 
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= j Vm^{X^) exp 



(A5) 



The Gaussian integration over $ in Eqs. ( |A2D , ( |A3| ) can be performed, resulting in 



Z{J} = const. Zb{J} / JJP^„(X„)e-^'=ff^*"^^ 



a=l 



where 



Zb{J} = exp 



(A6) 



(A7) 



with the bulk two-point correlation function Gh{r_,r^) corresponding to the action The 
effective action Ses is given by 

5eff{^, J} = \Y.j '^^-j ^-(^°) (A8) 

dX^J{r)G,{r,X^)^^{X^) . 

a Q 

Note that evaluation of Eq. ([A6|) requires functional integration over the curved manifolds 
Qa- This is facilitated by expressing the functional measure / T'\1'o(Xq) in terms of the local 
coordinates y, which itself comprise a flat manifold (the local coordinate system). To this 
end we introduce the new fields ipa{y) = ^^[-^^(y)]- However, this transformation requires 
some care regarding the integration measure dX^ in Eq. ( |A8|) as well as the functional 
measure / V^^{X^) in Eq. The resuh is [||] 



(A9) 

where the field 0a (y) = [ga{y)Y^'^'^a{Y) is given for each manifold Vta in terms of the 
determinant Qaij) of its induced metric 



The new effective action Ses is given by 

a/3 ^ 



(AlO) 



(All) 
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with the kernels 



A.p{y,y') = [ga{y)]'/'G,[X^{y),Xp{y')]My')Y/\ 
Wa{r,y) = G,[r,X^{y)][g^{y)f'. 



(A12a) 
(Al2b) 



The functional measure J T>(f)a (y ) on the right hand side of Eq. (|A9|) is the one conventionally 
used on a flat manifold. The corresponding Gaussian integrations can thus be performed, 
resulting in 



Z{J} = Zb{J} exp 
with the kernel 



(A13) 



^(z:,r') = J2jd''yJ ^V^.(z:,y)A;^(y,y')^/3(r',y') • (A14) 

Using A;^(y,y') = [(7a(y)]"'/'M,^(y, y') [f7/3(y')]"'/', where M,^(y,y') is the functional 
inverse of Gb[Xa{y), X^iy')] (with respect to both y, y' and the indices a, (3), one finds that 
the factors of [ga{y)Y'^ in Eq. (^) cancel. From Eqs. ([7^21) - ( [M4l ) one can thus read off 
the final result for the two-point correlation function, 

G(r,r') = G,(r,r') - / rf^^Z/ / G^r, X,(y)] M,^(y, y') G,[r', X^(y')] . (A15) 

Choosing = 1, corresponding to only one manifold, gives Eq. ( |2.3| ). 



APPENDIX B: SHORT DISTANCE EXPANSION OF THE STRESS TENSOR 

In this Appendix we consider the expansion of the two-point correlation function for 
a general massless field theory described by a Hamiltonian 7^{$}, to first order in the 
deformations of the height profile of a bounding surface. To this end, we introduce a new 
type of short- distance expansion of the stress tensor near a surface with the following scale- 
invariant boundary conditions: (a) the Dirichlet boundary condition $ = corresponding 
to the ordinary surface universality class, and (b) the boundary condition $ = oo describing 
critical adsorption, corresponding to the extraordinary universality class. 

A deformed surface 5* given by the height profile /i(x) (see Fig. P can be obtained from 
the fiat surface 5*0 with /i.(x) = by means of a coordinate transformation, which maps the 
space (x, z) on the space (x, ^) . We define this transformation by 
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X = X , z = z + /i(x) Q{z) , (Bl) 

where Q{z) is an arbitrary differentiable function with Q{z) = 1 for z < zq with some zq > 0, 
and which vanishes for z — > oo. We denote the Hamihonian with the flat surface 5*0 by 7i 
and the Hamiltonian with a deformed surface S by Ti.. According to the definition of the 



stress tensor Tj^ [^6|-^ the change of H generated by the coordinate transformation ( pi] 
can be written as 

d r 



n-n 



/„/^^E[^(Mx)e(--: 



T,fc(r) + 0{h'') , (B2) 



where HS denotes the half-space r = (x, z) with z > 0. Using the property ^j^. dkTu. = 
and the divergence theorem, one obtains 

n = n+ f d^xh{x)T,,{x,z = 0) + 0{h'^) . (B3) 

The contribution to first order in h is located at the (fiat) surface and does not depend on 
the specific choice of Q{z). The higher order contributions 0{h'^) cannot be transformed in 
this way, and will not be addressed in the following, ^^^(x, 0) = lim^^o ^zz(x, 5) represents 
a surface operator, which does not, however, need to be renormalized at the surface, so that 
its scaling dimension equals its canonical inverse length dimension of d [^ , pO| ] . 

In the following we consider the cumulant ($(r)$(r'))'-^ of the two-point correlation 
function in the system described by above the deformed surface S. Using Eq. (|B3|) 

one finds 

($(r)$(r'))^ = ($(r)$(r'))^ - j ci^x /i(x) (T,,(x, 0) $(r)$(r'))^ + 0{h^), (B4) 

where ( )f denotes the cumulant within the half-space HS bounded by the fiat surface 5*0. 
We now consider the limit p = \r\\ — rj|| — ^ oo (see Fig.|l]), so that we can use the short- 
distance expansion (SDE) of the order parameter $(r) near the surface. For the first term 
($(r)$(r'))^ in Eq. (|B4|) , the SDE is well-known: (a) for the Dirichlet boundary condition 
$ = 0, the SDE is given by |,| 

$(r||,2) = a;z(^l|-'^)/2_^<l>(r||,;z) +..., (B5a) 

OZ 2=0 

where 9^$(r||,z = 0) is a surface operator, r^n is a surface critical exponent, and a is a 
nonuniversal amplitude; (b) for the boundary condition $ = oo, the SDE has the form 
US 
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= ^ + feT^"T,,(r||,^ = 0) + ... , (B5b) 

where ($(r||,z))o is taken at the critical point of the field theory, and I is the identity 
operator. The amplitude Bt is universal. Equation ( [B5| ) in conjunction with the scaling 
behavior ($(r)$(r'))o' ~ p~^'^~'^~^^^ /(^, j) gives the result for a fiat surface 

($(r)$(r'))|^ ~ (^y)('?ii-'?)/2p-(rf-2+';ii) ^ p ^ oo . (B6) 

For boundary condition (b), one has r/y = d + 2 and the property (T^^) = has been 

used. 

For the second term J d^xh{x) (^^^(x, 0) $(r)$(r'))^ on the right hand side (rhs) of 



Eq. (p^), the above procedure cannot be applied directly because the integration of T^zi'^, 0) 
separates the points r and r' from the surface. To proceed, it is illustrative to consider first 
the case of a constant height field /i(x) = Hq. In this case, the integration of T3^(x, 0) simply 



amounts to a surface shift in the form ET 

ho I rf^x (r,,(x, 0) $(r)<|.(r'))o^ = + ^) " (^7) 

Consider for illustration the case for which only r = {r\\,z) is close to the surface, i.e., 
z <^ z' . Since ($(r)'^'(ZL'))o' small z behaves like a power in 2;, the z derivative on the rhs 
of Eq. ( |B7] ) is larger than the z' derivative by an amount of order z' j z. Now we recall that 
for the boundary conditions (a) and (b), correlations near the surface are suppressed, so that 
one can expect that on the left hand side of Eq. ( p7|) actually only a small integration region 
around ry contributes to the z derivative on the rhs. This suggests the operator product 
expansion 

T,,(x,0)$(r||,2;) = A(x - r,,, ^$(r|,, + ... (B8) 

for (x, 0) close to r = (ry, z), where A(x, 2;) is a representation of the delta function S^{x.) 
in D dimensions, i.e., 

rf^xA(x,2) = 1 , lim A(x,2) = 5^(x) . (B9) 

Note that dz^{r\\, z) on the rhs of Eq. (^^) is not a surface operator, since the z derivative 
is taken at a distance z > from the surface. The validity of Eq. ( [B8|) can be verified 
for various cases. For two-dimensional systems at criticality bounded by a line with the 
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boundary condition (a) or (b), it follows from the local form of the conformal Ward Identity 
For the Gaussian model with the boundary condition (a), it can easily be verified for 



any dimension d. For a model at criticality with boundary condition (b), Eq. ( P8|) is 
consistent with the form of (T2^(x, 0) $(r||, z))o known from conformal invariance arguments 
for any d with 2 < d < A [^,0]. For this system we have checked Eq. (|B^) also for the 



correlation function {^{r_)vil!))o with (p{r) = ^(l) ~ {^{l)) ^'^ ^^^^ {one loop) order in 
the $^ interaction. 

Let us go back to the second term J d^x h{x) {T^zi^, 0) $(r)'^(zi'))o' °f -^l- (|^1P 

with z and z' fixed. In order to obtain its leading behavior for p ^ oo, the x integration 
can be divided in two regions. Within one region, x is far away from both ry and rj|. Hence 
Eq. ( p5|) can be applied to both points ry and rj|. Within the complement region, x is either 
close to ry or to rj| so that Eq. (|B8| ) can be used. Due to the z derivative in Eq. ( pSf ) in 
conjunction with the scaling behavior quoted below Eq. (P3|), the contribution arising from 
the second integration region is by a factor p/z or p/z' larger than the contribution from the 
first integration region. Using Eq. (|B6|) , we conclude that the leading contribution for p ^ oo 
of the second term on the rhs of Eq. (|B4| ) is given by [A{r) + A{r^)]{^{r)^{r'))Q with the 
amplitude A{r) in Eq. ( |2.12| ). Thus we obtain the leading scaling behavior of {^{r)^{r'))'~^ 
for p oo quoted in Eq. 

APPENDIX C: STRUCTURE OF THE LOOP EXPANSION 

We consider the diagrams on the right hand side of Fig.^(b). According to Eq. (|5.2| ), 
in the (p, 5) representation, the distances 6q of the vertices [dots in Fig.0(b)] from the 
surface have to be integrated using d6o. To one loop order, only the three diagrams 
in the first line of Fig.0(b) exhibit short-distance singularities at 6q = 0. These diagrams 
consist of the following components: 

. = goip; 6, 5o) [see Eq. (2.8)] ; (CI) 

= e~P^° ; (C2) 

-LU = g,ip; 6, 6o) [see Eq. (4.4)] ; (C3) 

-Ll* = -p }C{p, 0) e-P^o + p }C{p, 6o) - Q e^f'^" (C4) 
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with the constant f2 = ^ /C(p, 5)s=o ; 



o 



A6'q-'' (C5) 



with the constant 



where the function Fi{6q) is regular for 6q —>■ 0. The constant B is given by 



B = D 



-1 



(27r: 



f/o(a)e-'" - 2f/(a)e-" , (C8) 



where f/o(p5) = U{p,0)/6 and f/(p5) = U{p,6)/6, with f/(p, 5) from Eq. (|;|). Note that 
in Eq. ( |C7| ) the terms in square brackets in Eq. ( [4.4| ), which correspond to the first hne in 
Eq. ( p.23|) , do not contribute. 

Writing A = Aq + eAi + 0{e'^) and B = Bq + eBi + 0{e^) with e = 4 - one finds 
that ^0 = ^0- This nontrivial fact is the reason why the 1/e poles due to the short- distance 
singularities of the first and the third diagram in the first line of Fig. ^ (b) cancel each other. 
The second diagram can be written as — Ce — In(jp)] » i i x + F2{p,6), with Euler's 

constant Ce and a pole- free function F2{p, 6). The l/e pole in this expression is then removed 
by the factor with Zi from Eq. (|5.8|) , that multiplies the zero loop contribution « i i x 



of the correlation function. The remaining, regular contributions, including those from the 
diagrams in the second line of Fig. 0(b), contain additional logarithmic terms in 6 which are 
not present if the surface was flat. One can then identify these logarithmic contributions, and 
show that they can be recast in the power law according to Eqs. ( |5.11 ) - ( ^.13 ). It should be 



noted, however, that here this exponentiation is not entirely based on an RG argument, but 
relies on the plausible assumption that the self-affine structure of the surface should result 
in pure power laws (without logarithmic corrections) for the decay of correlation functions. 

An analogous calculation leads to the quoted results for lateral correlations. Since in 
this case both points are located near the surface, only four of the six diagrams in Fig. ^(b) 
are different from each other. 
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